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Abstract. We present PHOENIX atmosphere models for metal-rich cool white 
dwarfs using improved line shapes for the Nal and KI resonance doublets. Pro- 
files for collisional broadening due to H2 and He based on the adiabatic represen- 
tation show strong deviations from Van der Waals interaction at short distances. 
Comparison with observed spectra that show extremely broadened Na I lines in- 
dicates that a He-rich atmospheric composition is required to explain the line 
strengths and spectral energy distributions. Our current synthetic spectra, using 
an expansion in powers of density to the third order optimised for brown dwarf 
atmosphere conditions, significantly underestimate the observed absorption in 
the far wings, even predicting smaller total line strength than a Lorentzian pro- 
file. This is shown to be due to the handling of multiple perturber interactions 
becoming inadequate for the extreme densities of the coolest white dwarfs. The 
density expansion would have to be extended at least to the 7'^ order for an 
accurate treatment of such conditions and might break down altogether in the 
densest objects. The results of a direct calculation of the unified profile should 
therefore be used for model atmospheres of cool metal-rich white dwarfs. Qual- 
itative comparison of the full adiabatic profile to the spectrum of WD2356— 209 
indicates good agreement with the observed line shape. Observations of the 
coolest white dwarfs may therefore serve as a laboratory for testing the physics 
of the deeper atmospheres and interiors of brown dwarfs and giant planets. 



1. Introduction 

Observations of th e oldest and coldest white dwarfs following the study of 



Oppenheimer et al.l ([200 L ) have revealed two stars showing very unusual wide 



and deep absorption at 5000 -6000 A. In WD2356— 209 this spectroscopic fea- 
ture is accompani ed by a unusual combination of blue B — V and extremely 
red V — I colours ( Salim et al.ll200j ). A second star with this absorption fea- 



ture has been identified in Data Release 2 of the Sloan Digital Sky Survey : 
SDSS J13 30 01.13-F64 35 23.8 (hereafter SDSS J1330-h6435 |HardIir2 llOOi), 
which also shows indications of an additional strong absorption around 4000 A. 
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These spectra bear intriguing resemblance to the strongly broadened Nal and 
K I resonance doublets observed in cool brown dwarfs of late L or T spectral type 
( Tsuji et al.lll999l : iBurrows et al.|[2000l ). Collisional broadening of the Nal lines 



in a metal-enriched, extremely dense white dwarf atmosphere thus lent itself as 
an explanation for these observations. 
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Figure 1 . Comparison of the observed spectrum of WD2356— 209 
(jOppenheimer et all 1200 1[ ) to PHOENIX spectra using different models of the 
Nal line profile. The fi r st one (smooth solid line) is calculated from the unified 
theory of lAllard etlHI (|1999D in its expansion to the third power in density, 
the other two are u sing a Lorentzian line shape with the half-widths given by 
lAUard et all (l2006h . 



In addition to the Nal line the spectrum of WD2356— 209 shows some ab- 
sorption features indicating the possible presence of other opacity sources. In 
a previous analysis we have attempted to model these feature s with molecu 



lar b ands expected to form in a cool hydrogen-rich environment (iHomeier et al, 



20051 ) . In particular, chemical equilibrium models predict MgH and CaH ab- 



sorption around 5000 and 7000 A, respectively, for effective temperatures below 
approximately 5000 K. However S/N and resolution of the spectra would not 
allow the identification of these bands with certainty. On the other hand strong 
arguments exist in favour of a helium-rich composition, the most important ones 
being the generally higher densities in the highly transparent cool helium atmo- 
spheres, which would easily provide the extreme perturber densities evident in 
the line profile, and much shorter diffusion timescales in hydrogen atmospheres 
( Dupuis et al.|[l993l ). making the presence of significant amounts of heavier ele- 
ments practically impossible without a continuous source of enrichment. 

Yet under these premises a helium-dominated atmosphere containing a 
smaller fraction of hydrog en-rich material coul d still provide a viable explanation 
of the observed spectra. iKoester et al.l ( 20051 ) have analysed an, albeit hotter, 
WD of such composition and found its combination of hydrogen lines with an 
unusually strong Balmer decrement and Ca, Mg and Fe lines best explained with 
the admixture of ca. 10~^'^ solar composition material to a helium atmosphere. 
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Such a model also obviates the need for removing the hy drogen from the ac- 
crete d material e. g. by means of the propeller mechanism (jWesemael &: Truran 
19821 ). 
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Figure 2. Number density of helium (left scale) and temperature profile 
(right scale) for 3000 K (lower curves) and 5000 K WD atmosphere models 
with varying H/He ratio 



2. Spectral Models 



Stella r atmosphere models were calculated using the PHOENIX code (jHauschildt Sz Baron 



19991 ). in its static, ID LTE setup. Formatio n of molecules in ful l chemical equi- 
librium is included by the equation of state of Allard et al.l ( 200ll ) as described in 
iHomeier et al.l ( 20051 ). though in the white dwarf models chemistry is dominated 



by only a few hydrogen compounds. This EOS does not take into account, on 
the other hand, non-ideal-gas effects relevant at extreme densities. But densities 
in the line formation region mostly do not exceed the limit of ~10^^cm~'^ (cf. 
Fig. [3]), where a fl uid-like state is appr oached and corrections for non- ideality 
must be included ( Kowalski et al.l [20051 ). 



Profiles of the alkali metal resonance lines for perturbations due to neutral 
He and H2 are calculated as the Fourier transform of the autocorr elation func- 
tion o f the d ipole moment within the adiabatic theory, as detailed in lAllard et al.l 
( 2OO3I . [2OO5I ) . In the PHOENIX implementation the line opacity is calculated by 



splitting the profile into a core component, which is describing the interactions 
at long distance with a Lorentzian line shape, and the far wings for close interac- 
tions that can produce detunings up to several 1000 A. The latter is computed in 
the low density limit using an expansion of the autocorrelation function in pow- 
ers of density, which for the present models has been developed to the third order 
and evaluated at a perturber density of lO^'^cm"^. This method, successfully 
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applied under the conditions of brown dwarf atmospheres, possibly compromises 
results at higher density ranges, where multiple perturber interactions become 
important even at close distances. 

2.1. Comparison with Observation 

Figure [1] shows the optical spectrum of WD2356— 209 in comparison to three 
of our models. In all cases shown here, the logarithmic ratio of He/H number 
densities is set to —2, i. e. ~ 1 % hydrogen with a corresponding amount of heavy 
elements at solar-composition ratios added to the helium atmosphere . 

A comparison of the WD2356— 209 spectrum with hydrogen/helium LTE 
atmospheres is presented in Fig. [TJ The strong dependence of the alkali lines on 
temperature, becoming quickly weaker with increasing Teg, is evident. For these 
models the spectrum becomes increasingly difficult to reproduce at temperatures 
above 4000 K, as the the width of the Nal doublet would require unrealistically 
high Na abundances, even up to the solar value. Such high metallicities would 
also conflict with the observed limits on molecular bands such as MgH and CaH 
and the K I doublet at 7700 A, all of which should be seen quite clearly if a solar 
abundance pattern is assumed. 
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Figure 3. Full unified profile of NaID2 perturbed by 10^^ cm^'^ of He at 
3000 K, compared to the 7*'' order expansion and Lorentzian 

2.2. Influence of the Line Profiles 

The synthetic spectra using line profiles calculated in the impact approxima- 
tion, thus including multiple-perturber interactions and yielding a Lorentzian 
lineshape, show that this simpler model can better reproduce the total strength 
and width of the Nal line, though the shape appears very poorly matched. One 
might therefore suspect that the expansion of the autocorrelation function has to 
be extended to higher powers of density for a correct implementation of the uni- 
fied profile. As illustrated in Figure [2] the perturber density in the line- forming 
region of an ultracool WD atmosphere can indeed exceed considerably the value 
of lO^'' cm~^ for which the present version of the far wing profiles was evaluated. 
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This is confirmed by directly comparing the hne strengths for higher den- 
sities. However Figure [3] shows that at a He density of 10^^ cm~^ even the 
7*^ order of the expansion falls short of the correct line strength in the wings 
at ^ 100 - 2000 A by almost an order of magnitude. This plot also illustrates 
how the unified prof ile approache s the i mpact approximation close to the line 
core as discussed bv lAllard et al.l ( 20061 ). At larger separations, i.e. for closer 
interactions, the deviations become severe, though the latter model at most 
wavelengths still matches more closely the correct line strength than the density 
expansion, especially in the red wing. At a few 100 A separation we find the 
third order expansion to underestimate the absorption by a factor of more than 
10, the 7**^ order still by about 4, but the Lorentzian profile by only 1 — 3. We 
can thus appreciate why the latter produces relatively better agreement with 
the observation. 




In Figure H] we follow the further transformation of the line profile when 
proceeding to yet higher densities. Ultimately a merging of the line core and 
the satellite to the blue of the original line position emerges, increasing the 
opacity around 5000 A by another factor of 10, while the red wing tends to 
asymptotically constant absorption which is actually fairly closely matched by 
the Lorentz profile (not shown). With these simulations we may now identify 
much, if not all, absorption shortward of the Nal line centre with the blue wing 
of the resonance line itself. Noting the sharp drop-off of opacity below 5000 A 
we can also expect that not much of the profile is formed at He densities above 
2 — 3 X 10^^ cm^^. It is tempting to trace the local flux minimum near 5000 A 
to the contribution of the satellite feature at higher layers, which is still well 
defined up to 10^^ cm~^, but an unambiguous identification should require both 
a full model atmosphere calculation incorporating these high density line shapes 
and observations at higher S/N. 
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3. Discussion 

We conclude that WD2356— 209 shows the effects of coUisional damping at signif- 
icantly higher perturber densities than encountered in the atmospheres of brown 
dwarfs. Detailed modelling of this object pushes broadening theory to new limits 
and thus provides a more stringent test of our calculations. This should advance 
our understanding of the physical conditions in ultracool, compact objects in 
general, and improve atmosphere models for brown dwarfs and other low-mass 
objects as well. Although these do not show such extreme broadening in the 
observable spectrum, the alkali lines provide significant opacity in the deeper 
layers of their atmospheres, thereby affecting the thermal structure at the con- 
vective boundary and possibly their cooling rate (I. Baraffe, priv. comm.). Study 
of ultracool white dwarfs like WD2356-209 or SDSS J1330-F6435 could thereby 
prove a test case for modelling the evolution of low-mass stars, brown dwarfs 
and giant planets. 
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